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Abstract 

An effect of tlie volume change upon proper ferroelastic (martensitic) phase 
transitions in cubic crystals is considered. Corresponding terms in the 
Ginzburg-Landau expansion of the Gibbs free energy are analyzed for the 
first- as weh as second-order phase transitions from cubic to tetragonal lat- 
tice under the action of uniaxial and hydrostatic pressure. The pressure effect 
on the critical temperature as well as on the phase transition anomalies of 
isothermal compressibility and linear thermal expansion coefficient are stud- 
ied and recent experimental data on thermal expansion anomalies in VaSi, 
In-Tl and Ni-Al are discussed. The non-linearity of thermal expansion leads 
to the special relation between the shear strain and volume change as a re- 
sult of the elastic energy minimization. This phenomenon can provide the 
transformation from FCC lattice to BCC one, observed in the iron alloys. 
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I. INTRODUCTION 



There is a growing attention in recent years to the physics of martensitic phase transfor- 
mations in metal alloys [0,0, though these phenomena were known by materials scientists for 
many years as diffusionless transformations characterized by specific transformation kinetics 
not described by the classical nucleation theory [Q. Two different kinds of martensitic 
transformations are known, i.e., an athermal and isothermal ones. In the athermal case 
the transformation begins at some start temperature Mg, but the parent phase still exists 
until the temperature goes down to Mf - a martensite finish point. In the isothermal case 
the transformation proceeds in time at a constant temperature and, generally speaking, 
could be completed in some finite time, which might be very long and depends in turn 
on temperature. There is a particular sub-class of a so-called "thermo-elastic" martensite 
within this class of transformation, which is characterized by the reversibility of a structure 
change - the alloy regains its high-temperature structure upon heating from a martensitic 
low-temperature state through the transition point. This means that the lattice of product 
phase is coherent with respect to the parent one. This phenomenon gives the ground for the 
"shape-memory effect" and is, thus, of a great practical importance. 

Due to a spontaneous strain release during the transformation some elastic distortions 
appear in the matrix of the parent phase surrounding the martensite nuclei. The mini- 
mization of elastic long-range energy determines the shape of new phase precipitates and 
provides an equilibrium state that has complicated heterogeneous structure which involves 
multiple-twin bands for deeper relaxation of the elastic strain ^ . Thus, a real development 
of transformation takes place in complicated conditions of apriory unknown external pres- 
sure from the lattice of the parent phase upon the regions under transformations. Thus, 
the heterogeneity of the system makes it difficult to analyze the thermodynamic of phase 
transition, and idealized single- crystal systems should be considered in order to study an 
equilibrium structure of the martensite phase as well as an equilibrium development of the 
t r ansf ormation . 
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The reversibility of the transformation means that a phase transition between equihbrium 
phases takes place and the structure coherence as well as the symmetry breaking at a critical 
temperature implies that the transition can be analyzed within the frame of the Landau 
theory. In this approach the difference in free energy between the parent and product 
phase is considered as a function of some order parameter which is equal to zero in a high- 
temperature symmetrical phase and becomes non-vanishing below the critical temperature 
0. The fundamental difference of the structure of a martensitic phase with respect to the 
parent lattice is well known to be spontaneous strain P], so this phase transition belongs to 
the class of proper ferroelastics p|. It means that the structure change takes place through 
the elastic instability of the crystalline lattice of a parent high-temperature phase with 
respect to a spontaneous homogeneous strain of a special kind 0. In other words, some 
combination of the elastic modulii vanishes at the critical temperature. 

It gives rise to the drastic lowering of the frequency of a corresponding mode of the 
acoustic vibration leading to a central peak of the inelastic neutron scattering [|10| and this 
is the reason why this sort of phase transition is often referring to as a "soft-mode" one. 
Though, the modes never become completely soften |]TT| and a finite (albeit - small) jump of 



the order parameter is usually observed, the martensitic phase transition still may be treated 
as a weakly- discontinuous and considered in the frame of the Landau theory of a continuous 
phase transition. A corresponding Ginzburg-Landau expansion of the free energy in series 
of the symmetrized strain components has been developed [p!^-[T^ and the heterogeneous 



fluctuations [|T5| , nucleation of the martensite phase around the defects |jT6|] and some other 



phenomena [jT^ were studied in such a formalism. 

In the present paper the role of hydrostatic strain in the martensitic phase transitions 
is considered for the case of a cubic lattice of a high-temperature phase. A cubic point 
symmetry leads to a coupling between the shear strain and volume change in the elastic 
energy expansion near the critical temperature, and the volume change can be considered 
as an additional order parameter not related with the symmetry breaking. The coupling is 
shown to lead to the transition anomaly in the thermal expansion coefficient as well as in 



the isothermal compressibihty. For the case of the first-order transition some effect is found 
to occur in the low-temperature phase even at the temperature region outside fluctuation- 
induced singularities. 

The effect of the uniaxial pressure, that preserves the symmetry of the low-temperature 
tetragonal phase is studied as well. It is an external field conjugated to the order-parameter 
which is known [0 to suppress the continuous (second-order) phase transition for arbitrary 
small field value, however, the weakly-discontinuous transition is preserved under the field 
lower than the critical one. The dependence of the transition temperature on the external 
uniaxial as well as hydrostatic pressure is derived and the critical pressure is found. 

The non-linear analysis of the volume change shows, that the elastic energy minimization 
can provide the FCC structure of low-temperature phase for the BCC parent lattice through 
the special relationship between the shear strain and volume change. Hence, non-linear elas- 
tic effects might be responsible for the challenging FCC - BCC martensitic transformation 
in pure iron and some ferrous alloys. 

The paper is organized as follows. A brief description of a weakly-first-order transition 
within the Landau theory is a content of the Section The Landau theory of a proper 
ferroelastic (martensitic) transformation is analyzed in Section |T^. The volume change due 
to such a phase transition is considered in Section Non-linear thermal expansion and its 
possible contribution to the martensitic phase transition is studied in Section 0. 

II. THE LANDAU THEORY OF THE FIRST-ORDER PHASE TRANSITION 

A. Continuous phase transition 

Let us recall briefly main ideas of the Landau theory of continuous phase transitions [0 . 
As Landau supposed, if the symmetry group of low-symmetry phase Gi is a subgroup of 
the symmetry group Gq of the high-symmetry one, than there is some variable t], called as 
an "order parameter" which is invariant under all the transformations from the Gi whereas 
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some transformations from Go change it. The thermodynamic potential - the Gibbs free 
energy could be expanded as a power series in t] near the critical temperature. As the 
thermodynamic potential should not change under the symmetry transformations which do 
not change the structure, the rj in the high-temperature phase should vanish. 

General expression for the (Ginzburg-Landau) expansion of the difference in Gibbs free 
energy between the phases has the following form 

Ag = ^{T-nW + ^v' (2.1) 

where is a critical temperature and the coefficients a and C should be positive. The 
equilibrium value of r] is determined by the minimization of AQ with respect to rj: 

dAg , d^Ag , , 

^— = and ^^>0. 2.2 

The solutions are the high-symmetry phase with rj = 0, stable for T > Tc and low- 
temperature phase with rf = a{Tc — T)/C that is stable for T < Tc. The dependence 
of on T is continuous in the critical point T^, hence this model describes the second-order 
phase transition. The Gibbs free energy as well as entropy changes continuously trough the 
transition temperature 

Ag = -^{T-T,f (2.3) 

but the derivative manifests discontinuity 

AC -T^^^ - r^'^Q _o^'T^ .2 5) 

- - - -20- ^^-^^ 



Vanishing of the coefficient of second degree in the Ginzburg-Landau expansion 
when the temperature approaches Tc leads to the critical fluctuations of the order parame- 
ter. Mean square of the homogeneous order parameter fluctuations is given by well-known 
expression 



{rf) oc 



(2.6) 



a|T-T,|-i 

Inhomogeneous fluctuations appear to be crucial in many cases, however, for the purposes 
of present study they can be neglected due to long-range nature of elastic interactions in 
solids PI . 



B. External field eff"ect 

Let us consider the effect of the external field E conjugated to the physical variable of the 
order parameter. In what follows it is an external pressure conjugated to the symmetrized 
spontaneous strain. The Ginzburg-Landau expansions takes the form 

A^ = |(T-T,y + ^r/^-Er/ (2.7) 

and the minimization of with respect to 77 leads to the cubic equation 

with discriminant 



a(T - T,)r] + Crf'-E = Q (2.8) 



(.0) 



3C ; V2C, 

It is seen that for any value of external field E the high-symmetry phase with 77 = no longer 
provides the stable solution of the Eq.( |2.2| ). Instead, we get 1] ^ for any temperature. It 
is known that the cubic equation has one solution in real numbers for Q > and three 
ones for the case of Q < 0. It means that for temperatures below some critical one that now 
depends on E 

T,{E)=T,~^(2CE^f (2.10) 

the additional minimum of the Gibbs free energy appears that corresponds to new phase. 
However, the initial phase described by the high-temperature solution of ( |2.8| ) provides the 
minimum with lower value of the free energy and is, thus, stable. The free energy behavior 
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as well as the order parameter dependencies on the temperature for different values of the 
external field are shown in the Fig.|l| and Fig.^. 

It is seen from the Fig.|l] and might be proven rigorously that different minima of the 
^Q{i]) curve have different energies for any temperature T < Tq. Thus, the high-temperature 
state remains stable throughout all the region of the phase co-existence. Only the condition 
of E = leads to the degeneracy with respect to sign of r] that implies the equal energies 
of different minima. It leads to a phase transition of the first order under the variation of 
external field at constant temperature T < Tc. In other words, the variation of temperature 
and external field act on the systems described by the expression ( |2.7|) in a different way, 
because the field variation may lead to the phase change but the temperature one may not. 

The external field suppress the transition, however some decrease in 

^(T, E) = a{T - T,) + 2Cr,\T, E) 

leads to enhanced fluctuations of the order parameter around Tc with smooth peak instead 
of divergence ( |2.(j| ) shown in the Fig.(|^). So, some transition anomalies are preserved in 
sufficiently small external fields E. 



C. Weakly-discontinuous phase transition in the Landau theory 

The first-order phase transition arises in the Landau theory when the symmetry of the 
system allows to have non- vanishing third-degree invariant composed by the order-parameter 
component |]7|,pO|. Corresponding term should be taken into account in the Ginzburg-Landau 



expansion: 

Ag = |(T - TcW + + ^V' - Ev. (2.11) 

Choosing the case of -B < 0, that gives positive t] in the low-temperature phase, we can 
write the free energy expansion in the following form 

A^ = — -A^ = -C'- — + — -aC, (2.12) 



with 



r 



52 ^ 



Tr) , and cr 



i. First-order transition in the absence of external field 



For the a = case minimization of the Gibbs free energy ( |2.12D with respect to ( imphes 
the equation 

^ = rC-C^ + C^ = 0. (2.13) 

For r > Tq = J there is the only minimum at C = corresponding to a high-temperature 
undistorted phase. The second minimum 

at C = i (1 + 

— 4r ) , or 

appears at To and corresponds to a low-temperature distorted phase which initially has 
higher free energy. The phase energies become equal at r.,, = |, though the supercooling of 
the high-temperature state as well as superheating of the low-temperature one are possible. 
It means that the first-order phase transition takes place at the temperature T^,. As well as 
in the second-order case, high-symmetry phase becomes unstable at r = 0. 

At the temperature of the first-order transition T,, the order parameter jumps from the 
?7 = to the 

2 B , , 

overcoming the activation energy barrier 

A5. = ^§;. (2,16) 
The entropy now has a finite change at transition temperature 



8 



that correspond to latent heat of the phase transition 

. ^ 2 

AQ = -^T.. (2.18) 

In order the weakly-first-order transition to be properly identified and clearly separated 
from the background of critical fluctuations around ( |2.ti| ), the transition jump of the order 
parameter should be greater than its mean fluctuation. In other word, the energy scale of the 
problem should be larger than the thermal fluctuation energy fc^T. It implies the condition 
for the value of the third-order coefficient in the Ginzburg-Landau expansion 

B > [c'kBT.f. (2.19) 

If this condition is not satisfied, the phase transition is "too weak" to be first-order and will 
be seen in the experiments as a continuous one. 



2. Effect of the external field on the weakly- discontinuous phase transition 



Substituting ^ = ^ + | into the Ginzburg-Landau expansion ( |2.12| ), we get the third-order 
term excluded [pl|] : 



Ag = ^C' + ^-ac + Ago, (2.20) 



where 



1 . r 2 

r = T : (T = (T 1 

3 3 27 

and A^o = — - —. (2.21) 

" 18 3 108 ^ ' 



This is equivalent to the free energy expansion (|2.71 ) for the second-order phase transition 
under the external field, the only difference consisting of the term AQq that is independent 
on (. It appears because the free energy is counted with respect to the {( = 0) or {( = — |) 
state, that implies A^o = AQ{( = 0) 7^ 0. 

The condition ( [2.2] ) leads to cubic equation ( ^.81 ) with the effective temperature f and 
field a instead of the real ones. The sign of discriminant 



Q = j + i^-j <xAa + 27a' - ISar - + 4r^ (2.22) 

of this equation indicates, whether it has one root or three ones in real numbers. The latter 
case corresponds to the appearance of different minima on AQ{(), second minima of the free 
energy appearing when Q(T,a) < 0. 

Hence, ( p^.20| ) can be considered as the Ginzburg-Landau expansion for the phase transi- 
tion between the states, related with different minima of the Gibbs free energy. The minima 
have non-zero values of order parameter, because the symmetry is broken already by the 
applied field for any temperature. As there is no symmetry breaking, it is not a true phase 
transition, described by the Landau theory, however, the undistorted phase with C = 
can be treated as an analog of ideal high-symmetry "praphase" [^] that would allow the 
symmetry reduction to both of the phases which provide minima of the free energy. It is 
interesting to note that the first-order phase transition in absence of external field appears 
to be equivalent to the second-order one under the action of 'effective' external field, the 
only feature of this situation is zero value of r] for one of two minima of /S.Q{Q. 

The phase diagram in (r, a) plane is shown at the Fig.^ Additional minimum of the free 
energy appears for cri < cr < cr2 with 

^i,2 = -^(l±(l-3r)i) + ^, (2.23) 



that leads to the hysteresis with respect to the external field Ao" = ^ (1 — 3 r 



According to an analogy with the second-order phase transition described by (|2.7|) , the 
different minima of the AG{C) have equal energy only at (T = 0. This is the condition of 
the first-order phase transition between corresponding phases and it determines the effect 
of applied field on the transition temperature 

n{a) = ?>a + \. (2.24) 

y 

For cr = we get naturally r*(0) = |. The Eq.( |2.24] ) corresponds to the straight line on 
(r, cr) plane (Fig.^). For r > | on this line the equilibrium phase has C = |- This state is an 
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analog of the undistorted high-symmetry phase of the Landau theory without an external 
field which is unstable for t < | and becomes a maximum of AQ, i.e. the energy barrier 
with a height of 

9 2 



4 6 324' 

for the first-order transition between two different minima with C,\^2 = ±V— 't, separated by 
the order parameter discontinuity 



AC = AC = -VI -27a . (2.25) 

As this line of the first-order transition in the phase diagram at (r, a) plane sep- 
arates states without symmetry-breaking relationship 0, it terminates in critical point 

= |,crc = ^). The discontinuity in order parameter as well as the potential barrier 
separating different minima of the free energy vanishes when approaching this critical point. 
There is no transition for a > cXc or t > r^, that means suppressing the weakly-first-order 
phase transition under the fields greater than the critical one. In a contrast with the second- 
order case where arbitrary small external field destroys the phase transition, here we find 
that the fields lower than preserve the transition. The critical point is an analog of 
the continuous phase transition from state with C = corresponding to the breaking of 
symmetry with respect to change of the C sign. 



III. PROPER FERROELASTIC PHASE TRANSITION 

A phase transition characterized by the appearance of spontaneous strain at the tran- 
sition temperature is called ferroelastic p[. When this spontaneous strain describes the 
symmetry breaking at the transition, and is, thus, an order parameter, the proper ferroelas- 
tic transition takes place. For the case of improper ferroelastics the spontaneous strain is a 
complimentary order parameter, coupled with the primary one. 

The free energy difference between the parent and product phases for the case of proper 
ferroelastic transition is due to the elastic strain and corresponds to the Ginzburg-Landau 
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expansion of the elastic energy in series of the strain components [0]. The second-order 
term in the Ginzburg-Landau expansion is a hnear combination of the second-order elastic 
constants that vanishes at the critical temperature. It is the eigenvalue of the lattice stiffness 
matrix corresponding to the relevant irreducible representation of the symmetry group of 
the high-temperature phase [[7|. The strain tensor components transforming with respect to 
this representation form the order parameter and the phase transition belongs to so-called 
"soft-mode" class, because it is accompanied by a noticeable softening of the corresponding 
acoustic mode of atomic vibrations [Q, visible as a central peak of the inelastic neutron 
scattering. 



A. Spontaneous strain in cubic lattice 



In what follows the case of cubic symmetry of a high-temperature phase will be consid- 



ered, that describes the martensitic transformations in the A15 compound |^3[ as well as in 
some metallic alloys. The spontaneous strain tensor has only diagonal components and the 



order parameter is composed by their symmetrized linear combinations [12 



V6 



^yy ~^ 



1 

V2' 



V2 = —[^XX - eyy) 



(3.1) 

(3.2) 



The rji corresponds to the extension along z axis without the volume change and 772 describes 
the strain non-tetragonality in XY plane. The critical acoustic mode is a transverse phonons 
distributing in (110) directions. 

The elastic free energy expansion can be written in the Ginzburg-Landau form 



(3.3) 



with the following combination of the elastic constants as the coefficients |T3| 



A a ,^ ^ 
2 2^ 



-(Cii - C12) 



(3.4) 
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Substituting 



B = ^(Ciii - 3^112 + 2C123) (3.5) 
C = — (Ciiii + 6C1112 — 3C1122 — 8C1123) • (3-6) 

4o 



rji = ?7sin^ and 772 = 77 cos 



we get the Ginzburg-Landau expansion in the form 



= |(r - T.) v'-jv' Mm + J v'- (3.7) 

The minimization with respect to 9 

-g^ = -Bv' cosm = and 

imphes sin(30) = ±1 depending on the sign of B. In the case of -B < we get three solutions 

(^,0), {-^7], ^7]) and {-^ri,-^ri) , (3.8) 

corresponding to the low-symmetry tetragonal phases obtained by the extension of parent 
cubic lattice along three coordinate axis. The free energy dependence on 77 for these minima is 
the single-component Ginzburg-Landau expansion ( p^.llD for the weakly-discontinuous phase 
transition. As the solutions are related through the symmetry transformation from the cubic 
point symmetry group and, thus, are completely equivalent, we can consider further only 
one of them, e.g. {rj, 0), without loss of generality. Hence, we can use the Landau theory for 
the case of single-component order parameter. 

The cubic symmetry allows the third-order term in the Ginzburg-Landau expansion to 
appear, hence, the Landau theory says that the phase transition should be of the first order. 
Indeed, for the case of Ni-Al and some other systems partial mode softening takes place and 
the finite strain appears at the transition, though the shear modulus decreases considerably 
near the phase transition temperature. Some other so-called "pre-transformation" phenom- 
ena were found in a number of alloys This case respects to the weakly first-order 



transition mentioned in Section |II C| . 
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However, this third-degree term appears to be very small in In-Tl, VsSi and some other 
alloys where the critical mode becomes almost complete soften, i.e. shear modulus 
(Cii — C12) vanishes as temperature goes to Tc. Central peak of the inelastic neutron 
scattering as well as other critical phenomena appear and the order parameter undergoes 



very small change at the critical temperature ||5D|. So, the transition is very weakly 
discontinuous, being sometimes of the second order within the experimental accuracy. The 
considerable enhancement in the acoustic wave magnitude occurs in the vicinity of the critical 
temperature Tc. Hence, this case corresponds to the the second-order phase transformation 



considered in Section II A. 



B. Effect of external stress on the transition 



An applied pressure gives rise the 'external' stress tensor E corresponding to linear term 



—eE in the free energy expression |^T|. The symmetry-breaking strain components are 
sensitive only to diagonal components of E, hence, relevant external stress is superposition 
of hydrostatic pressure P = Tt{E)/V3 with ones applied (uniaxially) along z axis 



El = i-E,, - Eyy + 2E,,)/V6 (3.9) 



and within XY plane 



1 

V2 



E2 — {ExX — Eyy). 



Hence, Ei and E2 are external fields, conjugated to the primary order parameter compo- 
nents Tji and 7^2, respectively, whereas hydrostatic pressure will be shown below affect the 
ferroelastic phase transition through the volume change 770 related with the order parameter 
by the coupling term proportional to r]Q{ril + rj^). 

Applied non- hydrostatic pressure Ei^2 breaks the symmetry of undistorted phase, moving 
out corresponding minimum of A^(?7i, 772) from the origin and lifts the degeneracy between 
three low-symmetry phases (|3.8|). The pressure which is co-aligned with the spontaneous 
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strain corresponding to one of these solutions, e.g. (?7i, 0), preserves the tetragonal symmetry 
of distorted phase and the pressure value Ei is an external field conjugated to the value of 
symmetrized strain rji as a single-component order parameter, considered in Section [11 C 2 



Otherwise the stable low-temperature state of the system has rhombohedral lattice with 
three different lattice parameters that is characterized by the pressure- dependent r]i and 
?72. In what follows the effect of uniaxial pressure Ei conjugated to t]i, which preserves the 
tetragonal symmetry of low-temperature phase is analyzed. 

In the agreement with the analysis of Section [II B| the uniaxial pressure was found ||32 



to suppress the ferroelastic phase transition from cubic to tetragonal lattice in the VsSi 
compound, where third-order term is very small and the transition is of the second-order. For 
Ni-Al alloy, where the first-order phase transition takes place, the uniaxial pressure appears 
p3| to shift the transition temperature linearly in complete agreement with the Eg. ( p.24[ ). 
It should be noted that though this alloy exhibit martensitic transition where spontaneous 



homogeneous strain is accompanied by so-called shuffle |Tl| related with q ^ critical mode, 
the central peak of the inelastic neutron scattering as well as noticeable softening of Cu—Cu 
elastic constant appear well above the transition temperature. Hence, this case can also be 
analysed in the frame of the Landau theory of ferroelastic phase transition. 



IV. VOLUME CHANGE IN THE ELASTIC ENERGY EXPANSION 

If the phase transition is sensitive to applied external hydrostatic pressure then there is 
a difference in the volume of elementary cell of the parent and product phases. The volume 
change, indeed, takes place in some cases and it was shown |^4[ that virtual volumetric 
strain could reduce the potential barrier for the system to overcome in the phase transition 
development. In order to analyze the volume change one needs to include corresponding 
terms into the Ginzburg-Landau expansion of the Gibbs free energy. Then the equilibrium 
state of the system should provide the minimum of the free energy with respect to both 
shear strain and volumetric one. 
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A. Linear energy of the thermal expansion 



In the Unear elasticity theory the trace of strain tensor is a measure of the volume 
change, the corresponding symmetrized linear combination of the strain components having 
a form 



Tr(e) 1 , 

VO = = ^{^xx + eyy + ezz)- (4.1) 



Due to thermal strain as well as external pressure, the expansion of the elastic energy in 



series of tiq should start with the linear terms |31 



miVo) = -noAo{T - TR)r]o + + PVo, (4.2) 

where Tr is some reference temperature for volume change, kq and Aq are the volume thermal 
expansion coefficient and bulk modulus, respectively and P is external pressure applied. The 
bulk modulus Aq = :^(C'ii + 2C12) is always positive. 

The minimization of A^o('7o) with respect to rjo implies equilibrium values of i]q{T) and 

m{T,p) 

r,o = Ko{T-To)-^ (4.3) 

A^o = - Y«o(r - Tof - ^ + KoP{T - To) (4.4) 

It is easy to see that usual thermodynamical expressions for isothermal compressibility 
(3t and thermal expansion coefficient 

_ d'AGo _ dvo_ I 

= dPdf = df ('-^^ 
are satisfied for such an expression of the Gibbs free energy if one bears in mind that P is 
opposite of the pressure inside the system, which is a thermodynamical variable (we also 
put the volume of whole system equal to 1 for convenience). 
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B. Effect of the volume change on ferroelastic phase transition 



The lowest-order term of coupling between the shear strains 771^2 and volume change rjo 



for the case of a cubic symmetry of the high-temperature phase is ^3 

AQrnt = D7]o{ril + rjl), (4.7) 
where D = ^^{Cm — C123). Substituting this term into the expansion of elastic energy in 



a power series of the volume change, we get the expression ( ^^ ) for the equihbrium value of 
?7o in the following form 

Vo = f^oiT -To)-^-^ivl + Vl) ■ (4.8) 

The coupling term does not induce any anisotropy in the {rii,rj2) plane, hence, three 
equivalent minima of the Gibbs free energy at low temperature have a form ( |3.8| ). We 
consider one particular solution of the form (771,0). The uniaxial pressure Ei applied along 
Z axis does not break the tetragonal symmetry of corresponding low-temperature phase. 
Thus, it is an external field conjugated to the rji order parameter and we can use the results 



of Section |1I C 2| to study its influence on the phase transition. The volume change r]Q as well 
as its derivatives - thermal expansion coefficient a and isothermal compressibility (3t does 
not depend directly on but only follows the dependence of symmetrized shear strain 771 
through Eqs.(|]8D. 

The general expression for the free energy has the following form 

/\g = AQo + Ag^ + AQint (4.9) 

where AQi is the Ginzburg-Landau expansion (|2.11|) with respect to r/i as the single- 



component order parameter with included effect of applied uniaxial pressure Ei (|3.9| ) 



A^i = f (T - T^Wi + f + - ^1^1 • (4.10) 



Substituting ( |4.8| ) with 772 = into Eq.(|4.9|) and taking Tc as a reference temperature Tr 
for the volume change, we get the following renormed Ginzburg-Landau expansion of AQ in 
the power series of rji 
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Ag{T, P, r/i) = AGo + ^l^i^ r^l + 
B' . C 



^vl + ^v'i-EiVi, (4.11) 



with the coefficients 



A'{T,P) = a' {T-T[) = 
2PD\ 

>i + 2koD){T - Te) - -^j , (4.12) 



B' = 5i and C' = Ci--— . (4.13) 

^0 



Here A^o {T, P) is the energy of a high-symmetry phase with ?7i = 0, given by the ( ^ . 
The critical temperature 

2DP 



is shifted by the apphed hydrostatic pressure 

dTr 2D 



(4.15) 



dP Ao{ai + 2KQDy 

in agreement with experimental studies, e.g. in In-Tl alloys |^^. New stiffness is given by 
the formula 

a' = ai + 2tiQD. (4.16) 

In order the Eq.( [4.11j ) be considered as an analog of Eq. (|2.11|) , the a' and C should be 
positive. 

The inequality C < Ci corresponds to the smoothening of the potential relief when 
additional degrees of freedom appear, which allows the system to relax easily. It should be 
noted that the inclusion of the terms corresponding to the possible volume change into the 
Ginzburg-Landau expansion ( [4.10| ) does not affect the third-order term hence, the order 



of phase transition could not be changed by applied hydrostatic pressure. 
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As the hydrostatic pressure change the transition temperature according to Eg. ( |4.15| ), 
the hne ( p.24| ) of the first-order phase transition becomes a surface in the 3D (T, P, Ei) phase 
diagram given by a formula 

nP,E,)=T, + ^ + ^P-^E,. (4.17) 

The transition can be induced by variation of hydrostatic pressure under the fixed values of 
temperature T and uniaxial pressure Ei at the point 

The critical point of the end of transition line ( |2.24| ) is now a line in the (T, P, Ei) phase 
diagram given by the uniaxial pressure E^. = — Pj'/(27C"^) and critical hydrostatic one that 
depends linearly on the temperature 

Ao ( Bl\ 



Pc{T) = a'{T - T, 



2P) V 3C' ) 

and vanishes when T goes to Tcp = Tc + Bf/{3a'C'). It means that some values of hydro- 
static pressure suppress the transition, caused by the change of uniaxial pressure Ei. The 
closer temperature is to the Tcp the smaller hydrostatic pressure needed for the transition 
to disappear. 

C. Transition anomalies of isothermal compressibility, thermal expansion coefficient 

and specific heat 

The isothermal compressibility, thermal expansion coefficient and specific heat are 
expressed by Eqs.f ^TSf ), ( [4.6|) and (|2.5|) , respectively, through the second derivatives of 
A^(T, P, ?7i) with respect to temperature and hydrostatic pressure for an equilibrium value 



of rji given by the condition ( p.2|) in the form 



J^{A', rii) =A'r]^ + Pi rjl + C ril - E, = (4.19) 



For the first derivative we have the following formula 
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A(A^_A^,).^^ + ^^ (4.20) 
dT ^ °^ OA' or dr]i dT ^ ' 

Second term in this expression vanishes as equihbrium r^i is given by the Eq.( |2.2| ) and taking 

into account Eg. ( [4.12D we get the formula 



^ (A^ _ A^o) = «' I , (4.21) 



that leads to the following expression for the transition anomaly of the specific heat 

ACp = -T^ ( - A^o) = T r^, (4.22) 
For isothermal compressibility and thermal expansion coefficient we get analogous expres- 



sions 



4^2 (9r)i 

A/^--^(A^-A^o)^-^..^ (4.23) 

A. = ^ (A^ - A^o) = r^i 1^ . (4.24) 

dPdT ^ °^ OA' ^ ' 



It is easy to see that the Keesom-Ehrenfest relationships 

dT, _ A^ _ TAa 
liP ~ 'A^ ~ ACp 

are satisfied. 

Differentiating both sides of Eq. (|4.19|) we get 

dr]i _ fdJ^y^ 
OA' ~ OA' \dr]i) 

and resolving ( |4.19| ) with respect to A', we can finally obtain the following expression 



(4.25) 



drii rjf^ 

dA'~ Ei + Bir]f + 2C'7]f 



* ^ = . „:S..^,„3 ■ (4.26) 



D. Second-order case 

Let us consider the case of the martensitic phase transition of the second order for which 
Bi = 0. In absence of the uniaxial pressure ( [4.2(j| ) does not depend on t]i and, hence, neither 
on temperature nor on hydrostatic pressure: 
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drji 1 



dA' 2C' ■ 

Low-temperature phase appears at with continuous evolution of the order parameter 
that gives the volume difference between the parent and product phases 

Ar/o = -^(T;-T) (4.27) 

and the discontinuities of the isothermal compressibility, thermal expansion coefficient and 
specific heat at phase transition take the following forms 

2 / D\'^ 2D^ 

C'Ao~ AoC,-2D^ ^^-^^^ 

- - ^'^ 2{A,C, - 2D^) ^^-^^^ 

Thus, to find three independent parameters of the model - a', C and we have four 
measurable values - ACp, A/3t and Aa along with 

^^^ = 1^ (4 31) 

dP a' 4) ^ ' 

which are related by Eq. ([4.25|) . 

As was mention above, for the case of second-order ferroelastic phase transition consid- 
erable critical fluctuations take place due to the softening of Cu — Cu shear modulus. The 
inhomogeneous fluctuations ?7i(x) appear to be relevant only in very close vicinity of the 
critical temperature |T^. From Eq.( ^.6| ) that describes the homogeneous order parameter 



fluctuations taking into account Eq.( |4.8| ) we get the following expressions for the critical 
fluctuations of volume 

77ooc-^|T-T,|-i (4.32) 
and thermal expansion coefficient 
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K= oc — T — Tc . 4.33) 

in the temperature interval around T^, where the fluctuations in rji are important. 

There are some experimental data available on the anomalies of thermal expansion in 
the single-crystal specimens near the martensitic phase transition which are in an 
agreement with the above results. In most of the alloys the thermal expansion coefficient 
increases near T^, that implies the positiveness of D. Outside the temperature region of 
thermal fluctuations the thermal expansion coefficient does not depend on the temperature. 
Besides, linear model of Eq. (^4.2| ) leads to the independence of a on the applied hydrostatic 
pressure. 

The second-order phase transition disappears under the applied external field Ei, hence, 
Cp, (3t and k manifest continuous behavior near the temperature under arbitrary small 
external field. However, from Eq. (^4.26 ) we get in such a case 

For sufficiently small values Ei of external uniaxial pressure the difference in the isothermal 
compressibility, thermal expansion coefficient and specific heat outside close vicinity of the 
transition temperature appears to be described by Eqs. ([4.28|) - (|4.30| ). The critical diver- 
gencies given by Eqs.( [4.32D and (f4.33|) disappear and the smeared peaks around appear 
instead. 



E. The first-order transition 



1. Absence of external uniaxial pressure 



If the third-order coefficient in the Ginzburg-Landau expansion (|4.1(]| ) has a non-zero 



value, then the Eq.( [4.1l[) in the absence of external uniaxial pressure describes the first- 



order phase transition at the temperature 



^* - + 9 a'C " + Aoa' + 9 a'C ^^'^^^ 
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The shift of the transition temperature by the apphed hydrostatic pressure has the same 
form ( f4.14| ) as for the second-order transition. The volume difference between the phases 



appears to depend on the temperature through the temperature dependence of rji given by 
Eq-CT 

1 \ 2 



This leads to finite volume change at the phase transition temperature 



^VoiTi) = (4.37) 



that can be observed in diffraction as well as dilatometric studies. 
Eq.( |4.26| ) takes the form 

dr]i r]i 



dA' Bi + 2C'r]i 



^ l+fl-^(^-^c')l 1 (4.38) 



2C' r ' V* Bf 

and along with Eq. (|4.24|) gives the temperature dependence of difference in thermal expan 
sion coefficients between the low- and high-symmetry phases in the form 



a' D ( Aa'C , „ \ - i , , 



The change of thermal expansion coefficient at the transition point T'^ now has a form 

The value of Ak decreases to that given by Eq. ( |4.29| ) as T goes down from T'^ to T'^. 

For the difference in isothermal compressibility between low-temperature phase and high- 
temperature one we can find similarly 

that gives us the phase transition discontinuity in the form 

23 



A/3t(T:) = 1^ . (4.42) 



It should always be positive according to general thermodynamical arguments 
The specific heat has the following temperature dependence 



Ta'^ ( ( Aa'C 



with the jump at the transition temperature 

ACp(r;) = 2T;^ = 2r;:M^^J^ (4.44) 

The phase transition discontinuities of volume and entropy are related through the 
Clapeyron-Clausius relationship with the slope ( |4.15| of the equilibrium line at phase di- 
agram 

dP AS a'Ao' ^ ■ ^ 

Both the thermal expansion coefficient and isothermal compressibility of undistorted 
phase with rj = do not depend on temperature, so, the expressions ( [4.39| ) and ( [4.41| ) de- 



scribe the temperature dependence of these quantities in the low-symmetry phase that can 
be observed experimentally below the transition temperature T^. It should be noted, how- 
ever, that the effects can be seen for being sufficiently far from T^, outside the temperature 



region where thermal fluctuations p.6| are important, because the fluctuation- induced singu- 
larities of the thermal expansion coefficient as well as the other quantities become larger in 
critical region than the jumps in their equilibrium values. Thus, the third-order coefficient 
Bi should satisfy condition (|2.19|) . 

Such a situation occurs in the case of Ni-Al and some other alloys where shear modulus 
at Tl is soften only slightly - by 10 -^ 20%, and the experimentally measured tempera- 



ture dependence of shear modulus can be interpreted as pointing even to negative JTl 
However, for In-Tl system where third-order term appears to be very small and the shear 
modulus almost vanishes at the transition temperature, the volume discontinuity is so small 



that it is hidden by thermal fluctuations ( [4.321 ). The similar effect occurs with respect 



to other discontinuities at the transition temperature T^, which is very close to Tc. 
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2. The effect of external uniaxial pressure on the anomalies around the first- order phase transition 



For the case of the first-order phase transition the dependence of the discontinuity in the 
order parameter on the uniaxial pressure Ei follows from general expression ( p^.25| ) 

that leads to the following volume change at 

A,„(Ti) = -- A(,;) ) (4.46) 

Both parent and product phases have //i 7^ and 7^ under applied external uniaxial 
pressure, thus, the differences in the isothermal compressibility, thermal expansion coefficient 
and specific heat are proportional to 

/ dr]i\ dr]i^i dr]i^2 



dA' J '^'^ OA' dA' ' 

where rji i{T, P) and ?7i,2(^; P) correspond to two different minima of the free energy given 
by the different solutions of Eq. ( [4.19[ ) . Using the dimensionless variables ( |2.21| ) we can write 
Eq.( [4.26D in the form 



dvi 1 (C + -'' 



r/i — ^ 



dA' C a-{C + i)2 + 2(C + 
At the transition point we have a = and 



1 



that gives the expression 



,,|!?L=K±i)!. (4.47) 

dA' 2C'C 

Taking into account that has the same value — f at the transition point for both phases, 
we can finally obtain 

A(r/i||)=-|7(l-27a)-. (4.48) 
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From the Eqs. ( [4.22| ) - ([4.24] ) we get the expressions for the phase transition discontinuities 



of the specific heat, isothermal compressibihty and thermal expansion coefficient as follows 



ACp(T:) = 2r:^fl + ^I^') ' (4.49) 



In the limit of small values of external uniaxial pressure we get the Eqs. ( [4.4^ ), (f4.42|) 



and ( [4.40|) , derived from their temperature dependence in the Section P^V E 1| . When Ei goes 



to the value of the critical point Ec, these discontinuities diverge as oc \Ec — E- 



V. TRANSFORMATION FROM FCC INTO BCC LATTICE VIA 

SPONTANEOUS STRAIN 

A. FCC BCC transformation through the Bain strain 

There is the case of martensitic transformation of especial interest, namely FCC - BCC 
transformation in Fe and some ferrous alloys. Since, there is no group-subgroup relationship 
for the symmetry breaking, the Landau theory is, generally speaking, inapplicable to this 
case. However, there is an orientational relationship between lattices of the parent and 
product phase, and the transition could be described in terms of spontaneous strain of 
so-called Bain type 0. 

The Bain strain is the single-axis shear of the same kind as an order parameter of 
the ferroelastic phase transition from cubic to tetragonal lattice. It is accompanied by the 
volume change, that is approximately 1.5% in the case of pure Fe, where transformation from 
austenite FCC 7-phase to ferrite BCC a-one takes place at 910°C. If the lattice periods 
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for austenite and martensite (ferrite) are and Oq, respectively, then the strain tensor 
components have the form 

exx = (^yy = ^f^— - ^ and e^z = — - ^ ■ (5.1) 

The fundamental feature of this compared with the above considered phase 

transition from cubic lattice to the tetragonal one, is the fixed value of the spontaneous 
strain needed to get the symmetry properties of the low-temperature phase. In the above 
considered case for any non-zero value of the order parameter t]i the symmetry of the lattice 
was tetragonal, whereas in the case of the FCC - BCC transformation the peculiar value of 
rji is needed to get the low-temperature BCC lattice. If we separate the shear strain from 
the volume change by taking the latter equal to zero, then we get single (and very large) 
value of symmetrized shear strain ( |3.1|) 



V6 V2-1 
"■ = - 272 + 1 "'"'^^^ 

Hence, this phase transition is completely different from the continuous ones, which the 
Landau theory describes, where the value of the order parameter changes with temperature 
in low-symmetry phase according to the minimization of its Gibbs free energy (|2.2|) . 

However, coupling with the volume change rjo makes it possible to have the rji variation 
in low-temperature phase without breaking of its symmetry. Indeed, the strain tensor ( p.l| ) 
implies the following expressions for the symmetrized combinations used above as the order 
parameter components 



(2v^2 + l) _^ 



Vo = ' ^ ' — -V3 (5.3) 



73 



a. 



7 



r^, = -JhV2-l)^ ; r/2 = 0, (5.4) 
V 3 

and we get the relationship between shear strain and volume change in the form 

2^2 + 1 ^ 
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Hence, the variation in value of rji preserves the BCC structure of the low-temperature 
phase, if rjo is changed in such a way that this relationship is satisfied. It should be noted 
that ( |5.5| ) is meaningful only in restricted region of rjo and rji. For example, r/i = implies 
unreal result Oq, = from ( |5.4| ). Thus, ( [5 .51 ) is justified only in some vicinity of the transition 
that is characterized by small volume change ^/3r]Q. 

Having supposed FCC - BCC transformation to be ferroelastic one, we should get the 
minimum of elastic energy for the values of rjo and t]i obeying the condition (|5.5|) . Let us 
study what are the coefficient in the expansion which provide such a minimum. Without 
careful analysis, it should be noted, however, that the linear approximation used above 
gives another kind of relationship ([4.8|) between the shear strain and volume change, thus, 
non-linear approximation for the thermal expansion energy should be used. 



B. Non-linear elasticity for the volume change 

Non-linearity arises naturally when taking into account large value of the strain tensor 
component. The Tr(e) for Bain strain in pure iron is approximately three times larger than 
real value of the volume change for this transformation given by direct multiplication of the 
lattice periods of low-temperature phase 

— = {l + e^^){l + eyy){l + e,,) -1. (5.6) 

It could be expressed through the symmetrized combinations 7]q, rji and 772 as follows 



Vo + Vo 

(5.7) 



2^3 3^6 

For the FCC - BCC transition we have 772 = and proper account for the volume 
change should, thus, involve the terms of higher order in 770 and rji. Terms of the first 
and second order in the volume change 770 within non-linear approximation have no longer 
simple relation with the thermal expansion and isothermal compressibility that was obtained 
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in Section [IV A| . Similarly, other terms in both AQi and AQint should be changed. General 



non-linear elastic energy expansion near the elastic instability with respect to rji now has 
the form 

Ag = Lor]o + Y Vo + Y ^0 + D Vo Vi 

where the coefficients Bq and Ci in highest-order terms should be positive and we again 
consider the particular expression with 772 = 0. 
The minimization with respect to rjo implies 

^ = Lo + AoVo + Bo Vo^ + Dr],^ = (5.9) 
that leads to the following relationship between rjo and 771 in distorted phase with r/i 7^ : 

For the high-symmetry phase we have 

As there is no co-existing high-symmetry equilibrium states with different values of rjo, the 
condition 

Al-ALoBo = (5.12) 
should be satisfied. Substituting this expression into Eq. (|5.10| ), we get 

In order the right-hand side of this expression to be positive condition D < must be 
satisfied, because of positiveness of Bq. 

Finally, we can get the following expression for the minimum of the free energy 
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and Eq.( |5.5| ) along with (|5.12| ) lead to the following relations between the coefficients in the 
elastic energy expansion ( |5.8| ) 

Lo = 3Bo (5.15) 
Ao = 2VSB0 (5.16) 
X,_^B„. (5.17) 

3-2V2 ^ ^ 

These relations could be, generally speaking, satisfied only in isolated points on phase dia- 
gram and the phenomenological approach used in the present study is unable to find their 
origin. It could be done only in some microscopic theory beyond the scope of the paper. 
However, as far as these relations are satisfied, we can try to find their consequences for the 
elastic properties of the system under phase transition. 

Substituting these expressions into the free energy (|5.8|) and excluding the volume change 
through Eq.( |5.5D we get renormed expansion of the elastic energy with respect to sym- 
metrized strain rji 



V 116-41V2 ; 3 4 ' ^ ^ 

which can be considered as a Ginzburg-Landau expansion for ferroelastic phase transition. 
The first term does not depend on rji, the critical temperature Tc is defined by the condition 

A + V3Bo'-^±^ = 0. 
3-2V2 

and the temperature of the first-order transition with finite jump in rji and r]Q is given 
by an equation 

9A^ -3AB^ + CiB^ = 

where A and B are the expressions in brackets of the first- and second-degree terms in 
Eq.( p.l8D . In order some transition line to exist on phase diagram, the Bq coefficient should 



be temperature- and pressure-dependent. 
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VI. CONCLUSIONS 



We have analyzed the volume change effect on ferroelastic (martensitic) phase transitions 
and considered the case of cubic lattice of high-symmetry phase as an example. The mini- 
mization of elastic energy with respect to hydrostatic strain as a secondary order parameter 
is shown to renorm the second- and fourth-order coefficients of the Ginzburg-Landau ex- 
pansion of elastic free energy in powers of symmetrized shear strain. The coupling between 
shear strain and volume change appears to shift the transition temperature under applied 
external hydrostatic pressure and lead to the finite volume effect of the weakly discontinuous 
ferroelastic phase transition. 

The isothermal compressibility as well as thermal expansion coefficient is shown to di- 
verge near the critical temperature of the second-order ferroelastic phase transition due to 
the homogeneous fiuctuations of the order parameter. The difference between their values in 
parent and product phases outside the fiuctuation region appears to be proportional to cou- 
pling coefficient. For the case of first-order transition isothermal compressibility and thermal 
expansion coefficient depend on the temperature in the low-symmetry phase according to 
the square root law. 

The uniaxial pressure conjugated to the symmetrized shear strain is shown to suppress 
the second-order transition, leading to the change of divergencies for smeared peaks in the 
temperature dependencies of isothermal compressibility and thermal expansion coefficient 
around critical temperature. We have found the first-order transition surface at the phase 
diagram in coordinates of the temperature and hydrostatic as well as uniaxial pressure. 
This terminates at the line of critical point and the uniaxial pressure of magnitude lower 
than critical, shifts the transition temperature, but preserves the transition. The critical 
hydrostatic pressure that suppresses phase transition has linear temperature dependence. 
The order parameter discontinuity and the volume effect diverge at the critical line as well 
as the difference in isothermal compressibility and thermal expansion coefficient between the 
parent and product phases. 
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The coupling between the volume change and shear strain is shown to lead to the FCC 
- BCC martensitic transformation for some special relations between the coefficients in the 
free energy expansion. Though some fixed value of the Bain strain is needed to get the 
low-temperature BCC lattice, the volume change as a secondary order parameter makes it 
possible to have some temperature variation of the shear strain preserving the BCC lattice 
and changing its period only. The non-linear expression for the elastic energy of thermal 
expansion is shown to lead to proper relation between the shear strain and volume change 
for the minima of elastic energy. 
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FIGURES 




FIG. 1. The dependence of the Gibbs energy on the order parameter 77 under the apphed field 
for different temperatures n > T2 > tq > T3 in the case of the second-order phase transition. 
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FIG. 2. The order parameter dependence on the temperature in various fields for the case of 
the second-order phase transition. Dashed hue corresponds to the absence of external field, a = 0. 
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FIG. 3. Mean square of the homogeneous order parameter fluctuations around Tc (r = 0) under 
different external fields. 
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FIG. 4. The region of the phase coexistence. The dashed hne corresponds to points of the 

first-order phase transition. It terminates in the critical point (tc = 1/3, ac = 1/27). 
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